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Abstract. The Helfrich-Canham bending energy is identified with a non-linear 
sigma model for a unit vector. The identification, however, is dependent on one 
additional constraint: that the unit vector be constrained to lie orthogonal to the 
surface. The presence of this constraint adds a source to the divergence of the 
stress tensor for this vector so that it is not conserved. The stress tensor which is 
conserved is identified and its conservation shown to reproduce the correct shape 
equation. 



PACS numbers: 87.16.Dg, 46.70.Hg 
Consider an embedded two-dimensional surface in three-dimensional space. The 
positive definite geometrical invariant, 



is a measure of the energy associated with the bending of the surface. The embedding 
is described by three functions X = X 1 of two variables £ a , = 1,2,3; a, b, ... = 

1,2). The extrinsic curvature tensor K a b is then defined by 



where e Q = <9X/<9£ a = V a X are the two tangent vectors to the surface and n is the 
unit normal. Indices are raised with g ab , the inverse of the induced metric g a b = e a -e^. 
dA = y/gd 2 ^ is the area measure induced on the surface by X, with g =det g a b- 

For a two-dimensional surface F{X] is invariant under ambient space conformal 
transformations. Since its introduction by Willmore pQ, it has cropped up in a number 
of contexts. These range from elasticity theory [2], to cellular biophysics where it 
models the bending energy of phospholipid membranes and is known as the Helfrich- 
Canham Hamiltonian 0E1E], to its Lorentzian generalization describing the action 
defined on the two-dimensional world sheet of a relativistic string propagating in 
spacetime, proposed by Polyakov to model colour flux tubes in QCD 00 El- We will 
couch our discussion in terms of Euclidean surfaces. 

The vanishing of the first variation of F[X\ under an infinitesimal deformation of 
the embedding functions, X — > X + <5X, produces the shape equation |91ll0j 



where V 2 is the surface Laplacian, and K = g ab K a b is the mean extrinsic curvature. 
We remark that the shape equation is of fourth order in derivatives of the shape 




(1) 



Kab = e b ■ V a n = K ba , 



(2) 




(3) 
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functions. Let us also note that, a fact that will be used below, the shape equation 
can be written as a conservation law [111 ITS] , 

V a f Q =0, (4) 

where the stress tensor f a is 

f a = K (K ab - \Kg ab ^j e b - (V a K) n . (5) 

In this brief report, we would like to explore the fact that it is possible to rewrite 
the quadratic appearing in Q in the form 

K ab K ab = (V a n) ■ (V a n) , (6) 

an identity which follows immediately from the definition of K a t and use of the 
completeness of the basis {e a ,n}, g ab e % a e\ = 6^ — n % n 3 . This identity associates 
bending energy with non-vanishing gradients of the normal vector; it appears to map 
the bending energy into a non-linear sigma model living on the curved geometry of 
the surface (see e.g. Ch. 14 of [T3|): 

F a [n] = \J dA [(V Q n) ■ (V Q n) + A(n 2 - 1)] , (7) 

with A a Lagrange multiplier that enforces the constraint that n be a unit vector. The 
Eulcr-Lagrange equations that follow from the vanishing of the first variation of F a [n] 
under n — * n + 5n are 

-V 2 n + An = 0, (8) 

together with the constraint n 2 = 1. Using this constraint, the Lagrange multiplier is 
identified as A = — (V n) • (V a n), so that the Euler-Lagrange equations take the form 

-V 2 n-(V Q n)-(V a n)n = 0. (9) 

These equations are integrable |14|: this fact is one motivation for investigating possible 
links with the bending energy 

If n happens to be the unit normal to some surface, then the fact that V a n = 
K a bg bc e c , and the Codazzi-Mainardi identities V a K ab = V K, allow us to rewrite the 
Eulcr-Lagrange equations in a geometrical form as V a K = 0. Clearly, this is not even 
a distant relative of the shape equation <0. It is of second order in derivatives of n, 
and of third order in derivatives of X. Despite this fact, unfortunately, this erronous 
identification has been made frequently in the literature; not by Polyakov, though, 
who was well aware of the pitfalls of a hasty identification (see Sect. 10.4 of jjj). 
What F a [n] fails to do is to capture the fact that n is normal to the surface. If n may 
point in any direction, then -F^n] describes a Heisenberg ferromagnet on the surface. 
On the other hand, if we wish to describe the bending energy, the unit vector n must 
also satisfy the constraints 

n-e a = 0. (10) 

If the surface is fixed n is also; if n is allowed to vary, the surface must respond 
accordingly. This leads to an interesting reformulation of the bending energy. We 
introduce the novel functional, with additional constraints, 

F[n,e B) X] = JdA i(V„n) • (V a n) + ^ (n 2 - 1) 

+ A a (n-e a ) + f a -(e a -V Q X)] . (11) 
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The normalization n 2 = 1 and orthogonality constraints (|10f> are implemented using 
the Lagrange multipliers A and A a , respectively. The latter constraints may be thought 
of as a frustration of the non-linear sigma model. They would be simple enough to 
implement if the e a were any two fixed vector fields. The fact that the vectors e a 
are the tangent vectors to the surface, however, couples n not only to them, but also 
through them to the embedding functions X themselves; this connection is captured 
in the final set of constraints in 1)11(1 . This model is a special case of a general 
construction introduced by one of us in where the issue of implementing the 
necessary integrability conditions needed for the surface to exist is sidestepped. 

The Euler-Lagrange equations for n, e a and X, that follow from the vanishing 
of the first variation of F[n,e a ,X], together are completely equivalent to the shape 
equation J2J. Moreover, the Lagrange multipliers f a appearing in I|11J) coincide with 
the stress tensor © - this justifies the abuse of notation. How the shape equation 
comes out is rather interesting. First, the Euler-Lagrange equation for the embedding 
functions X provides the conservation law (gj, V a P = 0, since X appears only in 
the last constraint. The equations for n, and e a determine the form of the Lagrange 
multipliers f a : for n we have 

- V 2 n + An + A Q e a = 0, (12) 

and for e a , when the constraints are enforced, 

i a =T ab e b - A Q n, (13) 

where we identify the stress tensor associated with an unconstrained non-linear sigma 
model for n on the background intrinsic geometry of the surface 

T ab = (V a n) ■ (V b n) - i<? ab (V c n) ■ (V c n) . (14) 

If A a = 0, (|12J) reproduces the Euler-Lagrange equations for a non- linear sigma model 
ijHJl. The normal projection of (fT2")l gives, as before, A = — (V a n) • (V a n) , or A is 
(minus) the energy density; taking into account the constraint 1)10(1 . its tangental 
counterparts give 

A a = e Q • V 2 n = V fc (e a • V 6 n) = V b K a b = V a K . (15) 

Modulo the constraints (|10|) and the definition of K ab , we find that T ab can be 
rewritten as 

T ab _ K ac R b _ l g ab K cd Kcd = K f R ab _ }_ g ab R \ ^g) 

To obtain the expression on the right we have used the Gauss-Codazzi equation for a 
two-dimensional surface, 

X -Kg ah = KK ab ~ K ac K c b , (17) 

which expresses the Ricci scalar induced by X, equal to twice the Gaussian curvature, 
in terms of a quadratic in K ab . Therefore, inserting l|15|) for A a and (|16fl for T ab 
into the Euler-Lagrange equation for e a ifTBl , we reproduce the stress tensor (0 . We 
identify A a as the normal stress and T ab as the tangential stress in the membrane. 

To arrive at the shape equation, let us return now to the conservation law Q. 
Using (|13|) , we can re-express the conservation law in terms of tangential and normal 
components: 

V a A a + K ab T ab = , (18) 
V a T afc - \ a K ab = . (19) 



Helfrich-Canham bending energy as a constrained non-linear sigma model 



4 



Using (|15fl for A a and (|16|) for T ab , it is immediate to see that the first, normal, 
component coincides with the shape equation ©. The coupling of n to the tangents 
provides a source to the sigma model stress tensor: it is no longer conserved. The 
source is the Lagrange multiplier A a implementing the constraint (IS) . When A a is 
substituted into i|19|) . it is identically satisfied, a feature that has its origin in the 
reparametrization invariance of the model. We note that the Lagrange multiplier A 
does not appear in these equations. We also remark on the fact that we do not need 
to know how K a b itself responds to a deformation of the surface in this presentation. 

What is the status of the Euler-Lagrange equation for n? Substituting A and 
A a into (|12() we obtain a purely kinematical statement about embedded geometries 
identifying what the projections of V 2 n are with respect to the basis vectors, e a and 
n: 

V 2 n= {V a K)e a -K ab K ab n. (20) 

To see this, just take a divergence of the Weingarten equations V a n = K ab e b using 
the Gauss equations V a e;, = — A^n to express V Q e;, as a normal vector. 

There is a second independent quadratic invariant in the extrinsic curvature, 
involving the square of the mean curvature, K 2 . However, it does not lend itself to 
a simple expression of the form ©. On the other hand, the fully contracted Gauss- 
Codazzi equation, 

K = K 2 -K ab K a \ (21) 

identifies the difference between the two quadratics with the scalar curvature, defined 
intrinsically i.e. independent of the normal n. For a two-dimensional surface, 
the corresponding invariant is the Gauss-Bonnet topological invariant. In higher 
dimensions, the difference is still independent of n. As a result, both A and A a are 
unchanged with respect to the values given above. 

We note, in this context, that the winding number of n is 

Q = h I dAeabe ^n l V a n 1 V b n k , (22) 

where e ab and t % ^ k are respectively the Levi-Civita tensors on the surface and in space 
|15|. Modulo 1|21|) . Q is just the Gauss-Bonnet invariant of the surface. Thus once the 
constraint is implemented, the winding number is fixed by the topology of the surface. 

To summarize, it has been shown that the coupling of a sigma model to the 
geometry constraining the unit vector n to lie normal to the surface converts it into 
the Helfrich-Canham model — involving only the surface geometry. A lot is known 
about both models. One expects that this identification might allow results from one 
model to be imported into the other. In particular, one would expect this formulation 
of the Helfrich-Canham model to be potentially useful in statistical mechanics. The 
functional -F[X, n, e a ] defined by (|llfl is quadratic in n. This suggests that, in the 
evaluation of the partition function, n may be integrated out. This possible virtue of 
the formulation presented in this paper will be the subject of future work. 
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